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Abstract 

We use the obser vational Hubble parameter data (OHD), both the latest observational dataset 
(jStern et al.l (|2010aD . referred to as SJVKS) and the simulated datasets, to constrain Lemaitre-Tolman- 
Bondi (LTB) void models. The necessity of the consistency check on OHD itself in the LTB cosmology 
is stressed. Three voids are chosen as test models and are constrained using the Union2 dataset of 
SN la as well as OHD. Despite their different parametrization, the results from our test models show 
some indicating similarities, e.g., the best-fit voids obtained from OHD are all considerably broader 
than those from SN la. Due to the small size of the SJVKS dataset, the constraints are not conclusive. 
The constraining power of the future OHD observations are therefore investigated, through a Figure 
of Merit (FoM) analysis based on the Monte Carlo simulated data. We found that, in the case that 
the future OHD become more consistent with SN la, the results from the test models are almost 
unanimous: 1) as many as 32 OHD data points at the SJVKS-like uncertainty level are needed to 
give a higher FoM than the Union2 dataset of SN la; 2) precise observation helps reduce this required 
number; 3) increasing the survey depth does not always increase the FoM. On the other hand, if 
the future OHD and the Union2 dataset keep favor different voids, in a similar manner as they do 
at present, the la confidence regions obtained from the two probes should finally separate. We test 
this conjecture and found that, the minimum observational requirement (the size of the dataset, the 
uncertainty level and the survey depth) for this inconsistency to emerge depends strongly on the void 
model. 

Subject headings: cosmology miscellaneous; statistical — methods: numerical 



1. INTRODUCTION 

Since the direct su pernovae rneasurement of the de- 
celera tion parameter (|Riess et al.lll998l : iPerlmutter et al.l 
Il999| ). the last decade has seen a solidification of the 
standard cosmological model, which has about 75% dark 
energy driving the acceleration of a flat Friedmann- 
Robertson- Walker (FRW) universe. A very long ar- 
ray of observati ons, including cosmic microwave back- 
ground (CMB) (iKomatsu et al.ll201lD. b arvon acoustic 
oscillations (BAO) (iPercival et all I2010D. weak lensing 
(ISchrabback et al.l 1201(1) 7" etc., are consistent with this 
model. 

Yet a radically different scenario has attracted consid- 
erable attention recently as an alternative model choice 
to fit the observations. The idea is to suppose we are very 
near to the center of a large spherical void, which is it- 
self embedded in an otherwise homogeneous background 
(jCelerierl 120001 ). The dynamics of this kind of mod- 
els can be described in Lcmaitre-Tolman-Bondi (LTB) 
metric, hence the term LTB models. Constraints from 
supernovae, CMB, BAO, a nd Hn have shown that this 
model c an not be ruled out jGarcia-Bellido fc Ha-ugboUd 
(l200l: iBiswas et al.l (poTof ). but see iMoss et al.l (|2011l ) 
and lZibin et al.l (j2008l )). As the cross-check from the con- 
sistency of the results based on independent evidences is 
an effective, if not t he only, way to pin dow n the correct 
cosmological model (jPeebles &: Ratrall2003f ) , more obser- 
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vations are needed to better discriminate the standard 
FRW and the LTB scenarios. 

Observational Hubble parameter data (OHD) has been 
me asured through the aging of passively evol ving galax- 
ies (iJimenez fc Loebll2002t ISimon et al.ll2005D and BAO 



effect ICaztahaga et al J120090. Recmtiy, a third method 
has been proposed to o btain OHD from the gr avitational 
wave standard sirens (jNishizawa et al.l 120 lit ) . OHD is 
shown to be consistent wi t h the standard co smologi- 
cal model ()Lin et al.l [2009t iStern et al.l l2010aD . Vari- 
ous ot her FRW models h a ve als o been tested agains t 
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potential of future OHD observations in constraining 
FRW models is explored via simulation in lMa fc ZhansJ 
()201lD . O nly recently wer e this dataset used to test LTB 
models in lFebruarv et al.l (I2010D (hereafter FLSC). 

The current work deals with a similar purpose as did in 
FLSC, as we believe more discussions are needed, both 
on the subtleties when using OHD in LTB cosmology, 
and on its potential constraining power as the dataset 
enlarges in the future. The investigations on these two 
aspects, as well as the current constraints on various void 
models, constitute the main content of this paper. 

We briefly review the LTB cosmology and introduce 
several test models in Section [2l In Section [3] we present 
the physics and assumptions involved in obtaining OHD. 
This analysis then guides our choice in picking the proper 
observable in LTB cosmology correspondent to the OHD 
dataset, and signifies the necessity for the consistency 
check on OHD itself in the LTB cosmology. Section [4] il- 
lustrates how, in the LTB cosmology, the likelihood anal- 
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ysis should be done when OHD is used; the best-fit void 
parameters obtained in this section are used in Section [5] 
to compare between the LTB void models and the stan- 
dard model. The potential constraining power of the fu- 
ture OHD observations is explored in Section |6l through 
a Figure of Merit analysis; in the same section we also 
discuss whether the future OHD dataset can be used, 
when combined with the Union2 dataset, to exclude the 
void models. We conclude in Section [71 

2. LTB MODELS 

The LTB solution to the Einstein field equations de- 
scribes the dynamics of a spherically symmetric dust uni- 
verse. We take a brief review of the basic formulas in LTB 
cosmology, then we introduce the test LTB models used 
in our analysis. 

2.1. Basics of LTB Cosmology 
The background LTB line-element is written as: 



ds' 



1 - k{r) 



dr^ + A^{r,t)dVl\ (1) 



where ' denotes d/dr, and k{r) is an arbitrary function of 
r, playing the role of spatial curvature. The stress-energy 
tensor of the mass source is given by: 



A' 



-pir,mSl- (2) 

Note that the FRW metric can be recovered by impos- 
ing A{r,t) = a{t)r and k{r) — kr^ . The Einstein field 
equations read: 



k{r) k'{r) 



A2 AA 

\2 



where ' denotes d/dt, and 



A-" + 2AA + k{r) = 0, 



A^ 



(3) 
(4) 

(5) 
(6) 



are the expansion rates at the transverse and longitudinal 
direction, respectively. Integrating EqHJ we get 



Hi = 



M{r) k{r) 



A3 



^2 



(7) 



where M{r) is another arbitrary function. The evolu- 
tion of A(r, t) can be obtained by integratin g Eg 171 the 
stand ard parametric solutions are as follows ( Moss et al.l 

ImiD: 



. / N Mir) , 
A{r,t) =— - coshr?), 



t~tB{r) 



k{r) 
M{r) 
\~k{r)f/- 



(8) 



: (sinhry — if) , 



for /c < 0; 



, M{r) ,^ 
A{r,t) =^"7^(1 - cosTy), 



t-tB{r) 



k{r) 
M(r) 
'fc(r)3/2 



(9) 



{rj — sinry). 



for fc > 0; and 



Air,t)^{^-^)'\t-tsir)f'\ (10) 

for A; = 0. Here ts is another arbitrary function, referred 
to as the 'bang time', due to the singularity behavior at 

t = tB. 

Substituting EqlTlinto EqlSl gives 



M'{r) 
~A^ 



(11) 



Evaluating EqlTl at present time gives the interdepen- 
dence between Mir) and k{r), which, with a notation 
similar to that in FRW cosmology, can be parameterized 
as: 



M{r) =i/i„(r)a„(r)A3(r), 
k{r) =i/io(r)[a„,(r)-l]A2(r), 



(12) 



where a subscript denotes the current value [t ~ to) of 
a quantity. Substituting Eq ll2l and A{t = to^r) = A^ir) 
to EqlHl(or EqU EqUUl depending on the value of fc(0)), 
one can get the age of the universe as a function of r, 
given by 

to - tB{r) = ^ \ , ^(a»), (13) 



^f±o(r)' 



where 



/ x — \^x sin 



T{x) = { 2/3 



(^-1) 



3/2 



/\ — x—x sinh 



x = \ 



(14) 



In order to compare the theory with the observation, 
one needs to associate the coordinates with the red- 
shifts. This can be done by noticing the red shift equation 
(|Alnes fc Amarzguiouill2006t lEnavistll2008D 



dln(l + z) _ i'(r,t) 
~ ^ A!{r,ty 

and the equation of dt and dr on the light cone 



dt 
dr 



A'ir,t) 

v/i-^W 



(15) 



(16) 



As a Cauchy problem the differential equations can 
be solved numerically for any redshift, once the pro- 
file of rim and H±o are set. We do the equation inte- 
ration with a modified open-sou rce software easyLTB 
Garcia-Bellido fc Haugb"oM 120081 ). where Ao{r) = r is 
chosen to fix the gauge. The calculation process for the 
cases when fi^ 7^ is as following. 



1) Start from t ~ to, r = 0; 

2) Solve for 77 from the second equation of EqElor Eql9l 

3) Obtain the expression of A, A and A' from the first 
equation of EqElor EqlQl 77' and 77 in the result can be 
expressed as functions of A using the second equation; 

4) Solve A' from EqEl EqEl EqU EqUIl and EqUl 
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5) Integrate Eg IT51 and Eg fTBl over z to get the new t{z) 
and 

6) Iterate steps [2115] tih the desired z. 

The original easyLTB does not support the void models 
with < 0. We modify the code to allow for negative 
Ofc. For fZfc = 0, it is not valid anymore to use the para- 
metric forms. However, as both EqlSland Eql9l converge 
to EqllOl when |r2fe| — > 0, in practice one can always ap- 
proximate the correct answer with the parametric form 
by setting f2fe to a small nonzero value: following the cal- 
culation process above, letting fi^ = 10^^ would cause a 
10~^ error in 77 and lead in turn to an 10~^ error in A for 
each integration step; the error will accumulate during 
the integration but should not be larger than 1% for any 
sensible number of integration steps (normally smaller 
than 100000). 

One can tell the angular diameter distance directly 
from the metric form, Eq[Tl 



dA{z)^ A{T{z),t{z)), 



(17) 



hence the luminosity distance (the distance duality re- 
lation holds as long as one stays with the metric the- 
ory of gravity and docs not consider any exotic physics 
(lEthcrington 1933: Bassett & Kimz 2004;)), 



dL(z) = (l + z)M(r(^),t(z)) 



(18) 



The LTB dynamics discussed above consider only mat- 
ter, so they break down at high redshifts where ra- 
diation become important. This is one reason why a 
large family of LTB models are voi ds that are embed- 
ded in an F RW background uni verse (jBiswas et al.ll2010l : 
iZibin et aL , 2008.; Moss et al . 20TT1 ). Moreover, out of 
the consideration of meeting the prediction of the in- 
flation, an Einstein- de Sitter universe is often chosen 
as this background (jGarcia- Bellido fc HaueboUd [20081: 
[February et all [2010? ) . The LTB models so constructed 
are referred to as EdS voids. This typ e of voids a r e said 
to be rul ed out as a whole a ccording to lMoss et al.l ([201 1[ ) 
(see also [Zibin et al.l ([2008| )). as those voids which fit the 
CMB data will give to o low a J?o that will co ntradict 
the local observation. [Clarkson h Regij ([2011[ ) found, 
however, that there are ways out if one takes into ac- 
count the radiation when joining the voids to the back- 
ground FRW universe, for instance, by introducing more 
relativistic degrees of freedoms, or by considering varied 
baryon fraction and/or baryon-photon ratio. 

Our aim in this work is mainly in exploring the con- 
straining power of OHD data, compared with SN la data, 
both of which arc low rcdshift observations. Therefore 
we do not have to take into account the role played by 
the radiation. Besides, whether the EdS voids are feasi- 
ble or not, as test models they could in any case serve as 
a test bench of the constraining power of different obser- 
vational probes. Thirdly, as mentioned above, the EdS 
background is consistent with inflation, and is therefore a 
theoretically conservative choice. For these reasons, the 
test models we use in this work are all EdS voids. Also, 
we assume there are no isocurvature modes, i.e., varia- 
tion in the baryon fraction or in the baryon-photon ra- 
tio, correlated with the void, because these will enhance 
the inhomogencity in the structure formation time which 



might make OHD itself invalid, not to mention using it 
to constrain models (sec the discussions in Section [3]). 

2.2. Void Models 

The radial inhomogencity of LTB solutions leaves an 
almost totally arbitrary form of the radial profiles of ilm , 
H±o and ts- For flm, the much discussed void- like pro- 
files are favored. This is because of the requirement that 
the LTB model should have an expansion rate decreasing 
with the distance to the center, which in turn resemble 
obscrvationally an expansion rate increasing with time 
in FRW models. 

Besides, the gradients in the bang time, ts, c orre- 
sponds to a curre ntly non- vanishing decaying mode ([Silld 
[l97l[Zibh][200l . This would imply a very inhomoge- 
neous early universe, hence violate inflation. More im- 
portantly for the current work, this will lead to great 
inhomogeneities in the galaxy formation time and make 
the OHD dataset invalid. So is must be a constant, and 
we set it to be 0. H±o for the desired void models is then 
given by substituting t{r) = to into Eq ll3l 



where 



H±o{r) = HoTin^), 



Ho = l/^Oj 



(19) 



(20) 



and is given by EqlTil 

With the relation between H±o and il„i deter- 
mined by Eg ll91 we need only fix the specific pro- 
file of Qm.. The Constrained^ G BH (CGBH) model 
([Garcia-Bellido fc Haugboll^ [2008| ) is chosen to be one 
of our test models, 



Qmir) ^ 1 + (no - I) 



1 - tanh[(r - ro)/2Ar] 



1 -I- tanh(ro/2Ar) 



(21) 



where flo describes the density at the symmetric center, 
ro is the characteristic size of the void, and Ar describes 
the steepness of the void near the edge. Figure [1] illus- 
trates profiles of ilm and H±o in the CGBH model, with 
Ho, ro, and flo fixed at 74 kms"^Mpc~\ 0.05, 2Gpc, re- 
spectively. Different colors (red, blue, brown) stand for 
different values of Ar/ro (0.1, 0.3, 0.9, respectively). It 
can be seen that small Ar/ro corresponds to an Qm pro- 
file that is flat at the origin, i.e., fl'^\r=o = 0, while large 
Ar/ro more pointed ones (ri^|r=o 7^ 0)- Because of this 
feature, and also the consideration that Ar should not be 
greater than unity, we use dr = Ar/ro as a free param- 
eter instead of Ar. Note that despite the notation, ^l„i 
does not reflect the exact proflle of the matter density. 
At a given cosmic time t* , the density proflle pm (r, t* ) is 
determined through EqlJl 

The value of 17',^ at the origin turns out to be of great 
importance. Defining the curvature as K = k/r"^, one 
could tell from EqlT^l and EqllOl that the r-derivative of 
K, K' , is directly related to O^. Indeed, it is straight- 
forward to show that an LTB model with homogeneous 
bang time is smooth at the center, K'\r=o = Oj if and 
only if r2J„|r=o = and ■;^^\r=o is finite. This kind of 
voids are referred to here as the 'smooth voids'. Particu- 
larly, CGBH void is a smooth void only when dr — )■ 0. It 
has been argued that the smooth voids are more natural. 
More importantly, the luminosity distance-redshift rela- 
tion in the smooth voids differs qualitatively from that 
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flm and -ffxo as functions of r, 


with fixed i?Oi f^Oi 



and different dr. Small dr is seen to lead to profiles that are fiat 
near the origin, and large dr more pointed ones. 

of the FRW models, which could be used to differenti- 
ate these two scena rios with the future supcrnovae data 
(jClifton et alJ[2008h . Keeping this result in mind, we 
choose the rest of our test models to be always smooth. 
From the discussion above, this is to look for models 
with that is differentiable at the origin. The first 
such model we choose has the simplest Gaussian profile: 

f^™(r) = l + (Oo-l)cxp(^-^^, (22) 

with the parameter set as {fio, T^^^ Hq}. The second 
smooth model is obtained through a cubic spline inter- 
polation. The constraints we set, given three r values 
ri = < r2 < r-^, for determining the Vim, profile, are 
that 

(a) a„(0) = no; 

(b) r2J„|r=o = 0, i.e., a smooth void; 

(c) fi,„(r2) = f2r2; 

(d) a„(r3) = 1; 

(e) r2J„|r=r3 = 0, which together with (d) makes this 
model an EdS void. 

The parameter set for the Spline model is then {fio; '^2, 
r2r2, ^3, Hq}- Iu piacticc, to make sure r2 is no larger 
than fa we use r2/r3 as a free parameter instead of r2. 

Throughout the paper the CGBH model is treated as 
our main model. Most of the calculation details will be 
described only for this model, although the main results 
for all three models will always be shown and discussed 
together. 

3. OBSERVATIONAL HUBBLE PARAMETER FROM 
PASSIVELY EVOLVING GALAXIES 

In this section, we firstly give a brief review on how 
OHD is obtained and utilized in the FRW scenario, then 
we will figure out the complexities that arise when mov- 
ing to the LTB scenario. As mentioned in the introduc- 
tion, there are currently two sources of OHD data points, 
those from the passively evolving galaxies and those from 



TABLE 1 
The SJVKS dataset 



z 


H(z) 






0.1 


69 ±12 


0.17 


83 ±8 


0.27 


77 ±14 


0.4 


95 ± 17 


0.48 


97 ±62 


0.88 


90 ±40 


0.9 


117 ± 23 


1.3 


168 ± 17 


1.43 


177 ± 18 


1.53 


140 ± 14 


1.75 


202 ± 40 



the large scale structure (BAO). The latter method de- 
pends on the detailed evolution of perturbations which 
is not well understood in th e LTB cosmology, although 
progr esses have been made (jZibinl l2008t iClarkson et al.l 
|2009( ). Therefore, the OHD discussed in this work refer 
exclusively to the former. 

3.1. FRW scenario 
The Hubble parameter in the FRW models is given by 



a I + z dtct ' 

where a is the scale parameter, dtct is the variation of 
the cosmic time due to a small change in the redshift dz. 
One can measure H directly th rough the differential age 
method (jJimenez fc Loebl[200l . We briefly review this 
idea below. 
For any galaxy one has, 

Tca{z) = Tf + Tga{z), (24) 

which simply states that the cosmic age Tca at redshift 
z equates the summation of the time for the galaxy to 
form, Tf, and the age of this galaxy, Tqa- The latter 
can be determined spectroscopically. If we could find a 
group of galaxies that share a uniform formation time, 
i.e., Tp =Const., we would then get a handle of dtct by 
simply measuring the age difference of those galaxies: 

dtct{z) = dTcA = dTp + dTGA{z) = dTGA{z). (25) 

The passively evolving galaxies can be identified by 
figuring out at every redshift the oldest galaxies which 
define the 'red envelop'. In this case, one assumes the old- 
est galaxies formed at the same time, which is of course 
a natural assumption in an FRW universe (one may call 
it the galaxy-formation version of the cosmic Coperni- 
can principle). By observing 24 massive galaxy clusters, 
determining the red envelop and dif ferentiating t he ages 
of the galaxies on the red env elop. IStern et al. (2 010a[) 
add two new points to those of iSimon et al. (2005,). The 
resulted dataset (hereafter the SJVKS dataset) is listed 
in Tabled 

3.2. LTB scenario 

There are two rates of expansion in LTB models, de- 
fined respectively in Eql5] and Eql6l It turns out that 
the longitudinal expansion rate H\\ has the same form 
as Eq[23l hence it corresponds to the observed H{z) in 
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Table [TJ This can be seen by combining EqlB] and EglTSl 
which gives: 

H.=--^4^. (26) 

Ambiguities arise, however, from Eg 1241 when deter- 
mining dtct ■ We have seen in last subsection that a same 
formation time of the oldest galaxies is the basic assump- 
tion in obtaining OHD. However, in LTB models where 
the universe has a considerable background inhomogene- 
ity, this assumption becomes unreasonable and has to be 
dropped. So one question that must be answered before 
any data analysis is whether OHD can still be regarded 
as a valid dataset in the LTB scenario. 

In answering this question, firstly we noticed that, al- 
beit the overall uniform-formation-age assumption, the 
validity of a OHD datapoint requires a same formation 
time only inside the redshift bin where it is obtained 
(Eql23l-Eql25l). This is to say, OHD data points are lo- 
cally defined and obtained, and therefore can be correct 
even when the global density, hence the formation time 
of the oldest galaxies, at different redshifts varies much. 
Secondly, in models' side, a standard viewpoint (referred 
to as the Onion approximation) treats the LTB void uni- 
verse as a group of thin shells structured together and 
that inside each of these sp herical shells the m atter is 
homogeneously distributed (jBiswas eFaI]|2003). If the 
bin size for determining OHD is small enough, and the 
background density varies only slightly inside it, we could 
identify these redshift bins with the homogeneous shells 
above and apply the Onion approximation. Thirdly, for 
the SJVKS dataset, the size of each bin is between 0.1 
and 0.15, where the first limit is so chosen that the age 
evolution between the t wo bins is larger th an the error in 
the age determination ([Simon et al.ll2005() . As the pre- 
cision of the age determination improves, we expect an 
even smaller bin size. 

It is by noticing the three points above that we suppose 
that OHD could be still used in LTB models. In order to 
be sure about this, however, one needs the exact knowl- 
edge about the thickness of the shell given the size of a 
redshift bin, as well as the steepness of the density profile 
at the time the oldest galaxies formed. These informa- 
tion will not be learned until one obtains the constraints 
on the model parameters. So in practice one has to first 
assume the Onion approximation holds, and check the 
validity of this assumption after the data fitting is done. 
The fitting result could be accepted only if the Onion 
approximation turns out to be valid. We will do this 
consistency check at the end of the next section. 

4. CONSTRAINTS FROM THE CURRENT OHD DATASET 

4.1. Likelihood 

We adopt statistics to determine the most likely 
values, as well as the confidence intervals for the param- 
eter sets used in our test models. Below are the main 
formulas using the CGBH model as an example. For a 



TABLE 2 

Scanning priors of CGBH model parameters 



dataset of {Hf} with errors {cTi}, 



X 



is defined as: 



ro 


dr Ho 


no 


[Gpc] 


[km sec"-"- Mpc"-"-] 




0.2-9 


0.1-0.9 65-85 


0.03-0.50 



TABLE 3 
Best-fit CGBH model parameters 



CGBH 




dr 


Ho 






[Gpc] 




[km sec"-'^ Mpc^^ 




OHD 


7.56 


0.21 


74 


0.058 


SN 


2.53 


0.70 


/ 


0.078 


OHD+SN 


3.60 


0.85 


74 


0.049 



The likelihood function for the parameters given the 
dataset can be obtained via Bayes' theorem, 

/:(i^o,^^o,r•o,Ar|{i^f}) 

_ p{{HP}\Ho, ro, Ar)p{Ho, ^o, ro, Ar) 



Pi{HF}) 
^p{{HP}\Ho,no,ro,Ar) 

oc exp ( — - 



(28) 



E 



[H||(z,;ffo,f^o,ro,Ar)-if,:"]2 (27) 



where the second equation comes from the fact that no 
priors constraints are imposed on the the data and that a 
uniform prior is assumed for the parameters inside their 
respective scanning ranges. Table [5] lists the parame- 
ter range scanned in the CGBH model, for the two void 
models see Appendix [Bl 

4.2. Constraints from SJVKS Data 

We apply the analysis to the SJVKS dataset of 
OHD, with the underlying model being our three test 
models. In Figure [2] we show the resulted Itr and 2a 
marginalized likelihoods for the CGBH model. As a com- 
parison, we also plot the likelihood contours for the latest 
Union2 dataset of SN la. Like in the FRW case, i^o in 
LTB models is treated as a nuisance parameter for SN 
la data (see Appendix lAl for a discussion on this point). 
Table [3] lists the 1— D best-fit parameters of the CGBH 
model. Since the exact fitting result of a particular model 
is of only secondary importance, we put the other two 
models' results in Appendix [B] just for reference. 

Neither dataset gives satisfactory constraints. Never- 
theless, several instructive clues can be observed from 
Figure [2j 1)0HD prefers a void profile that is consider- 
ably broader than that from SN la; 2)Both SN la and 
OHD prefer a void that is almost depleted at the ori- 
gin; 3) Aside from the nuisance parameter Hq, the SN la 
Union2 dataset is generally more constraining than the 
SJVKS dataset of OHD. 

Note that the 1-D likelihood of dr from OHD is far 
from being Gaussian and the best-fit value is seen not so 
well-defined. This in turn weakens the first point of the 
conclusions above. However, wc find that the same three 
conclusions hold for the two smooth void models too, de- 
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CGBH, OHD 
CGBH. SH 
Gauss, OHD 
Gauss, SH 
Spline, OHD 
Spine, SH 




Fig. 2. — Constraints on CGBH model parameters from the SJVKS datasct of OHD and the Union2 dataset of SN la. The red and brown 
filled regions correspond respectively to la and 2cr confidence regions obtained from the SJVKS dataset of OHD. Overplotted are the la 
and 2a confidence regions obtained from the Union2 dataset of SN la, in solid and dotted lines, respectively. 

fit parameters given by the SJVKS dataset and show 
in Figure m the density profile at different cosmic times. 
The inhomogeneity is seen to become more significant 
and the void gets deeper as the universe evolves, because 
of a higher expansion rate inside the void. In order to 
see more clearly how the density varies inside a small 
redshift bin of which we are more concerned, we plot in 
Figurc[5]tlie relative density variation inside a Az = 0.15 
bin. More specifically, we calculate on a given time slice 
Ap{r(z))/ p(r{z)), where p{r(z)) is the density value at 
r, with r converted from z using the best-fit parameters. 
As stated in the last section, each OHD data point is ob- 
tained using galaxies that are in a redshift bin no larger 
than 0.15; so Az = 0.15 is a conservative choice for our 
purpose. At Tqa = O.SGyr this relative change of the 
background density are globally less than 4% (see Fig- 
ure [S]) . If the oldest galaxies were formed around this 
time, which is not unreasonable for structures to form, 
the Onion approximation is a reasonable one. On the 
other hand, if these galaxies formed too late, the influ- 
ence of the background density variation on the structure 
formation may become too great to ignore. 

To determine this formation time. Eg 1241 is used and 
from which we have Tp = Tca{0) - rGA(O). The age of 
the oldest local galaxies, Tga{0), can be approximated by 
taking the Y -intercept of the red cnvelojiQ on Figure 1 1 of 
iStern et al.l ()2010af ). We find however, the outcome de- 
pends strongly on the models used to fit the galaxy spec- 
tra, which could sometimes yield a Tg^(O) greater than 
14Gyr. This number would defy all our best-fit void mod- 
els, as well as the lat est ACDM model result, 13.79Gyr 
(jKomatsu et al.ir2011[ ). because Tp would then be nega- 
tive. So finally, we could not know for sure whether the 





r[Gpc] 



Fig. 3. — Curvature profiles of the best-fit models from the SJVKS 
dataset of OHD and the Union2 dataset of SN la. OHD prefers 
voids with curvature changing slow near the origin; while the best- 
fit voids obtained from SN la are more cuspy. 

spite their totally different parameterisation. This might 
be a sign that the differences are real and are model- 
insensitive. We are particularly interested in the first 
point: if the future data of OHD and SN la both hold 
their different preferences of the void profiles, the smaller 
confidence regions they give would finally detach with 
each other. This would then be a sign that the void 
model is incorrect and should be excluded. In Section |6] 
we will give a rough prediction on this based on simulated 
OHD data. 

4.3. Validity check of the Onion approximation 

Now that we have obtained the best-fit void models, 
we can use them to check if the Onion approximation 
we have made indeed holds, i.e., if the universe can be 
treated as homogeneous inside the redshift bins where 
the OHD data points are obtained. We take the best- 



^ As we don't expect a homogeneous age of all the first galaxies in 
the LTB models, the red envelop does not have the simple relation 
with the cosmic time as it does in the FRW models. But the ages 
of the galaxies on the red envelop can still serve as a lower bound 
of the ages of the universe at the corresponding redshifts. 
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Onion approximation is valid or not. Only qualitatively 
could we say that the oldest galaxies observed seem to 
be formed very early, so very likely the SJVKS dataset 
is valid. More precise observations are needed to deter- 
mine Tca{0), hence Tp. As we need to use the best-fit 
models to generate simulated data, we leave this as an 
open problem and assume from now on that the Onion 
approximation is justified. 
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Fig. 4. — Density profiles at diff'crcnt cosmic time. As the voids 
get deeper, the background density become more inhomogeneous. 
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Fig. 5. — Relative density variation inside a redshift bin of 0.15 
as a function of the redshift. Ap(2:) = p(z -\- 0.15) — piz), where 
p{z) = p{r{z)) is the density value at r on a given time slice, and r 
is converted from z using the best-fit models. The dotted vertical 
lines indicate the data points of the SJVKS dataset, so one can tell 
approximately the relative background density variation inside the 
bins where each data point is obtained. 



5. COMPARISON BETWEEN BEST-FIT VOIDS AND FLAT 

ACDM 

5.1. AICc Analysis 

The theoretical prediction from different models at 
their best-fit parameters, as well as the observational 
data of OHD and SN la, are plotted in FigH and FiefHl 
respectively. The current observation of these two probes 
are not able to exclude either. To determine which model 




CGBH 

Gaussian 

Spline 

- - ACDM 



0.0 



0.5 



Fig. 6. — Hubble parameter as a function of the redshift. The 
prediction from the best-fit CGBH model based on the SJVKS 
dataset is plotted in red solid line and that from the fiat ACDM 
model is plotted in blue dashed line. Also plotted are the SJVKS 
dataset. 

is preferred by the data, a useful tool is the Akaike In- 
formation Criterion (AIC. lAkaik^ (flQTl V The AIC for 
a model is 

AIC = 2k- 2\nCM, (29) 

where k is the number of the parameters, and Cm is the 
maximized likelihood function. The preferred model is 
the one with the minimum AIC value, so not only does 
AIC reward a better goodness of fit, it also punishes a 
larger number of parameters. 

For OHD however, the number of the data points, 
n, is comparable to k, and the probability of ovcrfit- 
ting becomes considerable. A second-order correct ion of 
AIC, AlCr, should be used in such cases f Sugiurj 
iHurvich fc Tsailll989l : [Burnham fc Andersonli200i ): 



AICc = AIC 



2k{2k + l) 



1 



(30) 



AICc punishes the extra parameters more harshly, and 
it converges to AIC when n becomes large. Note that 
only the difference between the AICc is meaningful, so it 
is convenient to omit the constant in IuCm, which is the 
same for the same dataset, and express AICc in x^- 

2k{k + l) 



AICc 



2k 



y2. - 
A mill 



-fc- 1' 



(31) 



We list the AICc values for our LTB test models and 
flat ACDM in Table S] We discuss this rcsuh again us- 
ing the CGBH model as an example. Although having 
an AICc calculated from SN data that is fairly close to 
that of the flat ACDM, the CGBH model is consider- 
ably less supported by OHD. This is not to say that 
the CGBH model fits OHD worse than flat ACDM. On 
the contrary, the minimum from OHD dataset are 
A;^i„(CGBH) ^ 6.97 and xLnC^CDM) = 7.66, respec- 
tively, i.e., the former fits the data better. Therefore, of 
the three terms in Eq l311 the third contributes the most 
to the difference in AICc- As this term will get smaller 
as the dataset enlarges, the performance of CGBH in the 
AICc test against OHD is expected to get better in the 
future. 
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TABLE 4 

AICc VALUE FOR DIFFERENT MODELS 





CGBH 


Gaussian 


Spline 


flat ACDM 


OHD 


21.7 


16.4 


29.1 


13.2 


SN 


537.9 


536.9 


539.5 


535.1 



5.2. Effective Quantities 

To better visualize the evolution of the observable LTB 
universe, it is usually convenient to construct some 'effec- 
tive' parameters in a way originally used in FRW cosmol- 
ogy. The evolution of these parameters basically tells us 
what we would observe if we assume an FRW universe, 
while the real universe is, say, an LTB void. The best-fit 
model mentioned below will be the one determined from 
OHD and SN la combined. 

The first such parameter we discuss in this subsection 
is the effective deceleration parameter, defined to be 



-f 



(32) 



which in FRW cosmology would be just the familiar de- 
celeration parameter q. In Figure [7] we show q'^^ de- 
rived from the best-fit void models, compared against 
the deceleration parameter calculated from the best-fit 
flat ACDM model. The void models are seen to mimic 
observationally an universe that experiences an acceler- 
ation period between z = [0, 1]. 




CGBH 

Gaussian 

Spline 

ACDM 



-1.5L^ 

0.0 



0.5 



Fig. 8. — The eS'ective Dark Energy equation of state as a function 
of the redshift. The red solid line corresponds to the best-fit CGBH 
model, the blue solid line the Gaussian void, the brown solid line 
the Spline void, and the blue dashed line the best-fit fiat ACDM. 



where 



H\\{zf 

Hi, (0)2 



(34) 



The evolution of this effective EoS is plotted in Figure 
[51 Contrary to FLSC, we did not find a —1 crossing of 
the effective EoS in any of our test voids, which might 
be partly due to the more recent datasets of both SN la 
and OHD we have used. 




0.0 



0.5 



Fig. 7. — The effective deceleration parameter as a function of 
the redshift. The red solid line corresponds to the best-fit CGBH 
model, the blue solid line the Gaussian void, the brown solid line 
the Spline void, and the black dashed line the best-fit fiat ACDM. 

The second effective quantity is the effective equation 
of state (EoS) of the dark energy. Before the existence 
of the dark energy is unambiguously verified, one should 
always bear in mind that the dynamical properties of the 
dark energy people try to figure out today may be really 
an effective effect of, say, a background inhomogeneity. 
This effective parameter can be constructed again from 

f dhiQ , , 

-f + - (33) 



w 



,cff 



3 (M{1 + z) 



6. FORECAST OF CONSTRAINTS FROM THE FUTURE 
OHD DATASET 

We have seen in the last section that l)tlie two- 
dimensional confidence regions of the parameters from 
OHD and SN la are overlapped (consistent), 2)the two 
observational probes prefer characteristically different 
void profiles. In this section, we explore the situation 
in the future with Monte Carlo simulated OHD data. 
The first subsection is devoted to explaining the essen- 
tial ingredients of our simulation. 

To facilitate the comparison between SN la and OHD, 
we introduce for each void model a characteristic param- 
eter plane, which we refer to as CPP for brevity. The 
constraining power of one observational probe is quanti- 
fied via the area of the confidence region on this specific 
parameter plane. As we will find later, this is directly 
related to the definition of the figure of merit. Take 
CGBH model for example, the two parameters ro and 
dr contains most of the information one would like to 
know about the void profile, so tq — dr plane is chosen 
as its CPP. The CPPs of the Gaussian void and Spline 
void are not so obvious. In fact for the Gaussian void 
model, a better way might be to compare the length of 
the 1-Dim parameter confidence interval, and the Spline 
model the volume of the 3-Dim (£7^27 ''2, r3) parameter 
confidence region. Here for the sake of consistency, we 
look for the best parameter pair and choose the CPP for 
the Gaussian and the Spline void models to be vq — ilp 
and r2 — ^r2, respectively. 

Based on the simulated data two questions are inves- 
tigated, both of which are related to the CPP chosen 
above. In the following we describe them in turn. Firstly, 
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it is possible that when more OHD data points are ob- 
tained, they become more consistent with the SN la data. 
On the GPP, this would mean that the la and 2a con- 
fidence regions obtained from OHD become more over- 
lapped with those from SN la. One would like to compare 
in this case the constraining power of the two probes. 
The first question (Question-I) is therefore raised, that 
when the area of the OHD's confidence regions on the 
CPP become smaller than that of SN la. The stan- 
dard way to do this type of comparison is via the Figure 
of Merit (FoM). We calculate FoMs for our simulated 
datasets in the second subsection. 

Considering the different preference in the void profiles 
of OHD and SN la, the second question (Question-H) one 
would like to know through OHD data simulation is how 
many OHD data points are needed for the la confidence 
region on the CPP totally detached from that of the SN 
la, if the future OHD data points hold the same void- 
preference as the current ones. This will be the content 
of the third subsection. 

6.1. Monte Carlo Simulation 

The fiducial function of H{z). This function sets 
the expectation value of the simulated OHD data. For 
the purpose of Question-H, where the future data points 
are assumed to have the same void preference as the cur- 
rent ones, a straightforward fiducial function of H{z) 
comes from the best-fit void models obtained from the 
SJVKS datasct of OHD. 

In dealing with Qucstion-I, on the other hand, where 
we assume the future OHD and SN la 'converge' in their 
best-fit void profiles, hence the model parameters, one 
natural choice is again to use the best-fit void models; 
only this time the model parameters are given by OHD 
and SN la combined. The problem is that the constraints 
from the current datasets are far from tight, as we have 
seen in Figure[51 the best-fit parameters are not Gaussian 
in many cases, hence not quite well-defined. The (flat 
ACDM) model, on the other hand, is consistent with and 
well-constrained by all current observations, and there- 
fore generalize the best knowledge we have about the 
observational universe; a fiducial model based on this 
model would seem to be a proper choice. However, flat 
ACDM, as do any other FRW models, differs charac- 
teristically from the void models. Using flat ACDM as 
the fiducial model, one risks introducing a systematic er- 
ror, which would in turn cause a shift and/ or a shape 
change in the parameter confidence region (|Kim et al.l 
[20M) . hence a bias in the figure of merit (see the next 
subsection). Noticing the merits and demerits of each of 
the two choices, we obtain the center values of the simu- 
lated datasct from a mixture of them (jRobert fc Casellal 
120041 ). More specifically, statistically half of the data- 
points are drawn from the flat ACDM, with ilm = 0.27, 
Hq = 73, which is consistent with the 7- year Wilkinson 
Micr owave Anisotropy Probe (WMAP) and Hq observa- 
tion (jKomatsu et al.ll2011l: iRiess et al.ll2011[ ) ; and m gen- 
erating the other half of the data, the void model with 
the combined best-fit parameters is used as the fiducial 
function. 

The discussion of the following two ingredients of 
Monte Carlo simulation should be the same when ap- 
plied to Question-I and Question-H. Where the fiducial 
function is used for illustration, we use the best-fit void 



model (combined). 

The uncertainty model. This model generalizes the 
statistical and systematic uncertainties about a future 
observation. There is, however, no such specifications for 
future OHD survey, t o the authors' know ledge. To work 
around this difficulty, iMa fc Zhang! (j2011[ ) examined the 
SJVKS dataset and proposed a simple phenomenological 
uncertainty model. Here, in order to gain more control 
on the simulations, we adopt an even simpler strategy. 
We introduce a percentage, S, to describe the maximum 
uncertainty for each data point, namely, for a given fidu- 
cial value Hi, the expectation of the uncertainty ai is 
taken randomly from the interval [-Hid, Hid]. This un- 
certainty model is surely over-simplified: not all errors in 
the observations are redshift-dependent. Nonetheless, as 
shown in Figure [9l a J = 25% error line does act as an 
envelop of the current observational uncertainties, except 
at z = 0.48, 0.88, where the da ta points have extraordi- 
nary errors (IStern et al.ll2010bD . 




0.0 



Fig. 9. — The HiS line with 5 = 25% forms an envelop of the 
errors of SJVKS dataset. Only two data points at z = 0.48, 0.88 
are seen to be beyond the envelop. 

Once ai is obtained, the mid-point of the observational 
data will then be generated from a Gaussian distribu- 
tion J\f{Hi,af). One such sample is shown, against the 
SJVKS datasct, in Figure [TOl 

Sample size and survey depth. The survey depth 
and the sample size arc another two important aspects of 
an observation that influence much the parameter deter- 
mination. We denote them as -Zmax and Ns , respectively. 
Therefore we have in total three control parameters to 
describe a possible future observation (or to generate one 
sample of data points): {S, Zmax, Ns}. In practice, S takes 
three values: 25% (current precision), 15% (according 
to a predi ction of OHD from Atacama Cosmology Tele- 
scope, see lSimoii et al.l (I2005D), 3% ( now can be achieved 
at z 0.42. see lCrawford et al.l ()2010l )): z^ax takes also 
three values: 1.5, 2, 2.5, and Ng takes four: 64, 32, 16, 8. 

6.2. Figure of Merit 

We deflne the Figure of Merit (FoM) to be the recipro- 
cal area enclosed by the 2a marginalized likelihood con- 
tour on the CPP of each of our void models. Higher FoM 
would then mean tighter, hence better constraint. The 
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Fig. 10. — One snapshot realization of the OHD dataset, with the 
uncertainty level S = 25%. The SJVKS dataset (diamond with red 
errorbar) are plotted for comparison. 



TABLE 5 

FOM OF DIFFERENT VOID MODELS 





CGBH 


Gaussian 


Spline 


OHD 


0.183 


0.107 


2.508 


SN 


0.272 


0.712 


3.707 



current values of FoM from OHD and SN la for different 
void models are shown in Table [5l the values in the OHD 
row are seen to be less than their counterparts in the SN 
row. Question-I could be translated then to the 'FoM 
language' as: when will the future OHD dataset give out 
a FoM larger than that of the SN la Union2 dataset? 

To calculate the FoMs of the simulated datasets with 
fixed Zmax and S, three steps are taken in turn: 1) gen- 
erating a sample of 64 OHD data points that are evenly 
spaced in the rcdshift range [0.1, Zmax]; 2) calculating 
the FoM based on this sample; 3) Extracting subsamples 
with Ng — 32, Ng = 16, Ng = 8 and calculating their 
FoMs in turn. Repeat these three steps for, say, 100 
times, and we get for each Ng an array of 100 FoMs, the 
median of which is taken as the typical FoM. Then we it- 
erate the whole process for different {zmax, ^} and finally 
get a 3-D matrix, with the element being the typical FoM 
at each {znie,^,S, Ng}. 

The result is shown in Figure [TTJ The size of the 
symbol is proportional to the value of the FoM, which, 
if greater than that of the Union2 dataset of SN la, is 
plotted in red, otherwise blue. 

Three void models arc seen to share several very similar 
features, both qualitatively and quantitatively. Firstly, 
as one may intuitively assume, the FoM increases with 
the sample size Ng. At S = 25%, an uncertainty level 
similar to SJVKS, a sample of 32 data points will surpass 
SN Union2 for all choices of Zmax- 

The second point is not surprising either, that lowering 
the uncertainty level S will increase the FoM. At S = 15%, 
a sample of 16 OHD data points are sufficient to surpass 
SN Union2 for all choices of Zmax- If we could achieve an 
overall S = 3% precision, iV^ = 8 is already good enough 
to give a FoM larger than that of SN Union2. 



The third point deserves a bit more discussion. One 
might expect that a deeper survey will tighten the con- 
straints. This is not always true, however, judging from 
the value of the FoMs. For example, when Ng = 16 
and S = 25%, the FoM of OHD with z^ax = 1-5, 2.0, 2.5 
are 0.209,0.203,0.234, respectively. That is, higher z^ax 
might lead to a lower FoM. This phenomena exist in all 
three void models. One possible reason for this to hap- 
pen is that, for a given A^^, especially when Ng is small, a 
higher Zmax leads to a looser redshift coverage, which in 
turn leads to more freedom for the theoretical curve to 
vary between the adjacent data points. The extra free- 
dom in the theoretical curve is realized of course through 
the extra freedom of the model parameters, hence a lower 
FoM. 

6.3. 1(7 confidence region overlap 

As discussed above, Question-II is motivated by the 
difference in the best-fit models obtained from the 
SJVKS dataset of OHD, and that from Union2 dataset 
of SN la. The logic is that if this difference holds for the 
future data, then one could finally tell the inconsistency 
from the model parameter confidence regions obtained 
from OHD and SN la, for instance, their la confidence 
regions on the CPP will separate. 

To simulate this, we generate new OHD data using the 
SJVKS-induced best-fit void models as the fiducial func- 
tion. We 'mark' at first those parameters encompassed 
by the SN Union2-induced la contour on the CPP. Then 
for each simulated OHD dataset, we determine its like- 
lihood and count the number of the marked parameters 
lying inside its la confidence region. No- For each sim- 
ulation control parameter combination, we generate 100 
samples, obtain 100 NqS and log the median value. This 
median number then reflects the typical overlapping con- 
dition for a given {A'^,, (5, Zmax} combination, i.e., if it 
equals zero then very likely one could see the separation 
of the OHD- and SN la-induced la confidence regions on 
the CPP. 

The results from different void models are not so ho- 
mogeneous as in the last subsection. The Gaussian void 
model is the easiest for the inconsistency to emerge: most 
of the control parameter combinations except those with 
Ng = 8,6 = 25% have a median No equal to 0. For 
CGBH model, only at the largest dataset size (TVs = 64), 
the highest redshift reach (zmax = 2.5 and the highest 
precision {5 = 3%) is the median value of No zero. In 
the case of the Spline void. No is nonzero for all the 
control parameter combinations, that is, from r2 — ^lr2 
plane at least SN la and the OHD are always consistent. 
Therefore, the prediction about when the OHD observa- 
tion is powerful enough that, when combined with SN 
Union2 dataset, can help rule out the void models, de- 
pends on the specific void model. This is our answer to 
Question-H. 

7. CONCLUSIONS 

We have discussed in this work the use of OHD, both 
the present and in the future, as a probe to constrain the 
LTB models. 

We have firstly showed, that of the two expansion rates 
involved in LTB models, the one in the radial direction 
corresponds to the observable in OHD, which arc ob- 
tained from the passively evolving galaxies. The latest 
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Fig. 11. — The Figure of Merit (FoM) in the control parameter space. In these right-handed coordinate systems, the x-axes are the error 
level (5, the y-axes 2max and the z-axes the size of the datascts Ns- Length of the plus symbol (+) are proportional to the values of the 
FoM. The symbol is plotted in red, and the line connected to it from below is solid, if the FoM value is greater than that of the Union2 
dataset of SN la; blue and dashed otherwise. 



OHD compilation, also the one discussed in this work, is 
the SJV KS dataset updated and compiled in lStern et al.l 
(|2010al ). A subtlety is then discussed on using OHD in 
the LTB scenario: on one hand, a synchronized oldest- 
galaxies forming time through out the universe is the 
basic assumption behind the OHD observation; on the 
other hand, LTB models have large scale background in- 
homogeneities, and there is no good reason to believe the 
first galaxies formed simultaneously all over the universe. 

The problem is severe because if the dataset itself is 
problematic, surely it can not be used to constrain any 
model. We have examined the method in obtaining OHD 
and found that only the homogeneity inside the redshift 
bin (< 0.15), where each data point is obtained, is de- 
manded for the data to be valid. This is not a unreason- 
able assumption for the LTB models. In the standard 
Onion approximation, the LTB models are approximated 
as being composed by different spherical layers, inside 
which the matter is treated as homogeneously distributed 
and the dynamics are FRW-like. However, whether a 
redshift bin can be viewed as one such layer should be 
checked given the model parameters. We argue there- 
fore a complete parameter constraining analysis should 
have a consistency check procedure after the standard 
analysis, and we illustrated this process on our test 



models with the SJVKS data at the end of Section |4l 
The result is inconclusive due to the large uncertainties 
in the absolute ages of the oldest galaxies. 

The constraints obtained from the likelihood analysis 
is indicating. SN la performs generically better than 
OHD judging from the likelihood confidence regions, but 
the two observational probes both favor the voids with 
the center largely depleted of matter. Besides, and more 
remarkable still, although the la confidence regions from 
OHD overlap with that from the SN la data in all of the 
test models, the former favors the broad and fiat voids, 
while the latter prefers voids that are smaller and more 
steep around the center. 

In Section [5] we compared the best-fit void models 
with the best-fit flat ACDM, firstly through the corrected 
Akaike Information Criterion (AICc) analysis. The lat- 
ter, according to the AICc value, is better than any of 
our void models. Indeed, especially for the complicated 
models, the small size of the SJVKS dataset of OHD 
greatly degrade their performance in the AICc analysis. 

The dynamical properties are then compared, in terms 
of the effective deceleration parameter {q'^^) and the ef- 
fective equation of state (EoS) of dark energy {w'^^). 
is seen to enter the negative region recently (z < 1) for 
all the void models, mimicking the recent acceleration of 
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the universe. of the void models are negative today, 
mimicking the dark energy's repulsive power. 

This general consistence and the characteristic differ- 
ence in the best-fit voids of SN la and OHD, obtained 
in Section m are taken further by using the Monte Carlo 
simulated data in Section[S] More specifically, we use the 
simulated data to predict when OHD become more con- 
straining than the Union2 dataset of the SN la (Question- 
I), and when the two probes become inconsistent with 
each other (Question-II), respectively. We make use of 
the Figure of Merit (FoM) analysis for Question-I, and 
the answer turned out to be surprisingly homogeneous 
for the three models: Ns — 32 is needed for the FoM of 
OHD to be greater than that of Union2 dataset, when 
S = 25%; N, = 16 when 6 = 15%; and A^^ = 8 when 
S = 3%; increasing Zmax does not necessarily increase 
the constraining power. 

The answer to Question-II is not so simple and is only 
semi-quantitative: to observe an inconsistency for the 
Gaussian void model, a dataset with slightly more data 



points than SJVKS {Ng = 16), and at the current error 
level {S = 25%) will do the job; but only when Ng, 6, z^ax 
are all at their most demanding level (64, 3%, 2.5, respec- 
tively) does the inconsistency show up; in the model with 
the most parameters, the Spline void model, the simu- 
lated OHD and SN la are consistent for all the possible 
{A^s, (5, Zmax} combinations. 
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APPENDIX 
A. Ho AS A NORMALIZATION FACTOR 



Gaussian 



Spline 






45 (normalized) 
90 (normalized) 



Fig. 12. — Normalized fi ~ z relation at different Hq. The distance modulus /i given by Hq = H are normalzed by multiplying the 
corresponding luminosity distances by H/H^^^ , where H^^'^ are the best-fit Hq value for each model, other parameters all take the 
best-fit values determined from the Union2 dataset. The difference between curves of H = 45, 75, 90 are well within the observational 
uncertainties. 

The observed distance modulus data is given by (IGuv et al.ll2007l : lAmanullah eraII[2nTnh 

= nig - M + ■ xi - /3 ■ c, (Al) 

where m^, xi, c are derived from the fit to the supernovae lightcurves and therefore are known before fitting the 
cosmology, M, ax and /3 are fitted together with the cosmological parameters to minimize the residuals in the Hubble 
diagram. 

The theoretical distance modulus /i are derived from the luminosity distance. 



In FRW cosmology, is proportional to the reciprocal of the Hubble parameter, 

d™^ = (H^^'^r^DUO), (A3) 

where Dl is a function of the cosmological components, 6. One can see how Hq^^ degenerates totally in this case 
with the absolute magnitudes of the supernovae in Ea lAll bv writing explicitly the statistic: 

_^ m*g~ 51ogio j^L(g, z,) - 25 + (-M + a ■ - /3 ■ c -f ^log^oH™"") 

^2 



Testing LTB with OHD 
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This makes it impossible to obtain constraints of from, say, the Union2 dataset in the FRW framework. 

In LTB models, the expansion rate goes into the distance in a more complicated way. Keeping ri(r) profile unchanged 
and varying 7?o, the pattern of distance modulus-redshifts curve also changes, as is argued and shown in Figure 4 of 
FLSC. However, we found that in the rcdshifts range we are concerned, z < 2.5, treating Hq as a normahzation factor 
in all our testing models leads to an error that is too small to be of any importance. 

More specifically, we calculate the luminosity distance for each model at their best-fit parameters, then redo the 
calculation with the value of Hq changed to be 45 (and then 90) and other parameters unchanged. The 'normalized' 
luminosity distance are calculated to be 

jN45 _j45 4^ 
"L — "L rrBST ' 

90 (^5) 

"L — "L J/-BST' 

respectively, where the su pers cript 'EST' denotes the best-fit value. The normalized distance modulus fi^aie then 
obtained according to Eq. IA2I The two specific number, 45 and 90, are so chosen that the values of Hq scanned in our 
calculation all lie between them, therefore this choice should serve as an upper bound of the possible error introduced 
by the approximation. The relative difference in the normalized distance modulus varying Hq from 45 to 90 is 

= ^N45 • 

The maximum value of (5^ at z < 2.5 for the CGBH model is 0.54%, Gaussian void 0.26%, and Spline void 0.05%, 
respectively, i.e., all less than 1% and can not be differentiated by the current observation (see Figure [T^] for an 
illustration). Therefore, we conclude that, at least for the models discussed and the rcdshifts concerned in the current 
work. Hq can still be regarded as a normalization factor. 

B. CONSTRAINTS ON THE GAUSSIAN AND THE SPLINE VOID MODEL 



2 




60 72 35 0.03 0.17 0.30 2 21 40 

H, iJo ^0 



Fig. 13. — Constraints on Gaussian void model parameters from OHD and SN la. The red and brown filled regions correspond respectively 
to 1(T and 2(7 confidence regions obtained from the SJVKS dataset of OHD. Overplotted are the Icr and 2(t confidence regions obtained 
from the Union2 dataset of SN la, in solid and dotted lines, respectively. 



The constraints on the Gaussian and the Spline void models are shown in Figure [T3] and Figure [HI respectively. The 
scanning ranges of the parameters are the same as and can be told from the axis in the plots. From the constraints 
one can tell why the answers to Question-II in Section [S] are so different for the Gaussian void and the Spline void 
models: in the Gaussian void model the la regions of vq — fio from OHD and SN la are only slightly overlapped, so 
a few more data points of OHD will be able to separate them; the confidence regions of r2 — ftr2 of the Spline void 
model, on the other hand, contain a much larger mutual part, therefore many more data points and/or much higher 
precision arc needed to separate them. 
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